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ABSTRACT 


The interaction between species has been a long standing scientific research both for mathematical ecologists and environmental 
statisticians, biologists or scientists to mention a few. In the present study activity we have indentified the full potential of using a 
numerical simulation which is computationally efficient to differentiate types of stability due to a variation of a co- existence steady 
state solution. The novel results which we have achieved on the implementation of this method have not been seen elsewhere; they 
are presented and discussed in this paper. 
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1. INTRODUCTION 


The interaction between two populations dates back to the earlier formulation of Lotka-Volterra system of continuous non- linear 


first order ordinary differential equations. In the theory of competing species, the co-existence steady state solution plays a 
significant role in the survival of two competing species (Burden and Faires, 2001). Since such a steady state solution is a point in the 
phase plane, varying one co-ordinate and fixing the other co-ordinate can have an impact on the type of stability for the two 
competing species. This proposed idea cannot be successfully tackled using an analytical method. Therefore, for the purpose of this 
study, we have utilized the method of a numerical simulation to quantify the impact of one co-ordinate of co-existence steady state 
solution on the type of stability. 


2. MATHEMATICAL FORMULATION 


Following (Murray, 2001, Yan and Ekaka-a, 2001, Ford et al, 2010, Kot, 2001, Beeby, 1993). We have considered the following 
system of continuous non-linear first order differential equation; 


ae = x—I1(t)(4—0.0003x — 1(r)—0.0004.x — 2(7) : 
a = x—2(t)(2—0.0002x —1(t)—0.0001x — 2(t) " 
t 


Subject to the ordinary initial conditions¥1(0) = Xy9 > O and x2(0) = x29 > 0 


Method of Analysis 

We have obtained a co-existence steady state solution (8000, 4000) by using the standard analytical method of Crammer’s Rule to 
solve the above equations (1) and (2). This steady state solution indicates that for the two competing species to survive together, x; 
species will have a population size of 8000 whereas the x2 species will have a population size of 4000. For the purpose of this study, 
we first fixed x2 population size and vary the x; population size and using this variations to study the stability of the expected co- 
existence steady state solution. Secondly, we fixed x; population size and vary x2 population size and similarly study the stability of 
the co-existence steady state solution. The results that we have obtained on the application of this method are presented and 
discussed below. 


3. RESULTS AND DISCUSSIONS 


The results that we have obtained upon the implementation of the above method of analysis are fully presented and discussed in 
this section. 


Table1 Evaluating the impact of *1 between 800 and 4400 on the type of stability 


Example |x 1 |x 2 | Eigenvalue A | Eigenvalue A2 | Type of stability 
1 800;}4000}2 150 5/0 809 5|/Unstabtle 
2 1200/4000/2 03 6 71]0 809 5|/Unstabtle 
3 1600/4000/1 9 2 8 4]0 3 916/Unstabtle 
4 2000|4000)]0 175 4} 1 8246]/Unstable 
5 2400 |4000 1 724 5)-0.0445)/Unstabtle 
6 2800/4000]1.6274)]-0.2674/Unstable 
7 3200}/4000]-0.4927}]1 53 27|)]Unstable 
8 3600}4000]-0.7200}]1 4 0 0/Unstabtle 
9 4000;4000/-0. 489} 1 3 4 9/Unstable 
1 0}4400/4000]-1. 792) 1 2 5 2/Unstable 


What can we learn from Table 1? From Table 1, the ten empirical examples of the two steady state solutions(&_1 x2) are 
unanimously unstable having two positive eigen values and eigen values of opposite signs. By the mathematical concept of a 
quantitative behavior of solution trajectories, two positive eigen values contribute to the unbound growth of the solution trajectories 
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(hence the steady state solution is said to be unstable). In this context, two eigen values of opposite signs signify that the positive 


eigen value is growing unbound and faster than the negative eigen value that is contributing to the decaying behavior of the 


solution trajectories. A similar observation has been made for Table 2. 


Table 2 Evaluating the impact of *1 between 4800 and 8400 on the type of stability 


However, 


Example | X 1 |X 2 | Eigenvalue Ai | Eigenvalue Az | Type of stability 
1 1/4800/4000;/-1.4106)1 170 6;]Unstabtle 
1 2/5200/4000]- 1.6429) 1 0829/Unstabtle 
1 3/5600/4000]-1.87611|0 9961/Unstabtle 
1 41/6000/4000)/;-2.1100/0 9100/Unstabtle 
1 51/6400/4000/-2.3444)]0 8244/Unstabtle 
1 6/6800/4000]-2.5794)]0 73 94/Unstabitle 
1 74/7200/4000]- 2.814 8)]0 654 8/Unstabtle 
1 8}/7600/4000]- 3.05060 570 6]/Unstable 
1 9/8000/4000]- 3.28 68]0 486 8]Unstabtle 
2 0;/8400/4000]- 3.523 31]0 403 3}/Unstabtle 


the co-existence steady state solutions continue to be unstable up to the steady state solution (10000, 4000), after 


which the instability is lost (Table 3 & 4). From our analysis, we have observed that the instability is lost between the steady state 
solutions (10000, 4000) and (10200, 4000). After this bifurcation phase, each steady state solution is consistently stable having two 


negative eigen values that contribute to the decaying behavior of the solution trajectories of x 


Table 3 Evaluating the impact of 1 between 8800 and 11800 on the type of stability 
Example | x 1 |X 2 | Eigenvalue Ai | Eigenvalue Az2 | Type of stability 
2 1/8800/4000]- 3.7600]0 3 2 0 0}Unstabdtle 
2 2/9200/4000;/-3.99 7 8]0 23 7 OJUnstabdtle 
2 3/9600]/4000;- 4.23 42)0 15 42;)Unstabdtle 
2 4/10000/4000]/-4.4716]0 .0 7 1 6/Unstabtle 
2 5}10200}/4000|}-4.7091]- 0 0109);S t a b |! _e 
2 6|/10400/4000/;- 4.946 8]- 0 093 2;S t a b |! _e 
2 7/10800/4000;-5.1847/)]- 0 175 3}/S t ab! _e 
2 8/11200]/4000;-5.4227]- 0 25 73/;/;S t a bil e 
2 9/11400]/4000);- 5.660 8]- 0 3392;);S t a b |! _e 
3 0/;/11800)/4000;-5.8991]-0.4209);S t ab le 


(t) and x2(t). 


Table 4 Evaluating the impact of 1 between 12800 and 16400 on the type of stability 


Example |x 1 |X 2 | Eigenvalue Ai | Eigenvalue Az2 | Type of stability 
3 1|12800/4000]- 6 13 74]-0.5 02 6;/S t ab lee 
3 2/13200|]4000)- 6 3-7 5° 8 90 5 842);S t a ble 
3 3/13600|]4000)- 6 6 14 3/]- 0 665 7/S t ab le 
3 4}14000/4000]- 6 852 9/- 0 7471/S table 
3 5|14400]4000]- 7 0 9 15 )= 0 828 5/S t ab le 
3 6/14800|]4000/- 7 3 3 0 2/- 0 909 8/S t ab iloe 
3 7/15200|/4000)- 7 5 69 0);- 0 9910/S t abitloe 
3 8/15600|]/4000)- 7 807 8/- 1 0722/S t able 
3 9;/16000|]4000)- 8 0467/- 1 153 3/S t ab |! _e 
4 0;16400|]4000)- 8 2.8 5 7 |) 2 4 23 43/S t able 
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What is our next task? Having varied x; and fixed x2 to construct several steady state solutions that were unstable (dominant) and 


stable in the minor, we would like to explore further on the variation of x2 while x; is fixed. The results of this new idea are shown in 
Tables 5-8. 


Table 5 Evaluating the impact of *2 between 400 and 2200 on the type of stability 


Example | x 1 |X 2 | Eigenvalue A; | Eigenvalue Az | Type of stability 
1 8000/4 0 0);-1.13581]0 495 8/Unstable 
2 8000/6 0 0);-1.2853]0 5 25 3/Unstable 
3 8000/8 0 0);-1.42711]0 5 47 1/Unstable 
4 8000/1000);-1.5630]0 5 63 O0}/Unstable 
5 8000/1200);-1.69421]0 5 74 2/Unstable 
6 8000/1400);-1.8215]0 5 8 15/Unstable 
7 8000/1600);-1.9455]0 5 85 5/Unstable 
8 8000/1800);-2.0668]0 5 8 5 6/Unstable 
9 8000/2000);-2.1856]0 5 8 5 6/Unstable 
1 0;/8000/2200)-2.3024)]0 5 8 24/Unstable 


Table 6 Evaluating impact of *2 between 2400 and 4200 On the type of stability 


Example | x 1 |x 2 | Eigenvalue 41 | Eigenvalue Az | Type of stability 
1 1/8000/2400/-2.41721]0 5 77 2}]Unstable 
1 2/8000)/2600);)-2.5304)]0 5 70 4}/Unstable 
1 3}/8000)/2800)-2.6420]0 5 62 0]Unstable 
1 4/8000/3000/-2.75231]0 5 52 3}/Unstable 
1 5}/8000/3200/-2.8613]0 5 413}]Unstable 
1 6/8000)/3400]-2.96921]0 5 29 2}/Unstable 
1 7/8000)/3600);)-3.07601]0 5 76 0}Unstable 
1 8}8000)/3800]-3.1818]0 5 01 8;]Unstable 
1 9}/8000);4000);)-3.2868]0 486 8/Unstable 
2 0};8000);4200);)-3.3909]0 470 9/]Unstable 
Table 7 Evaluating the impact of *2 between 4400 and 6200 on the type of stability 
Example | x 41) X 2 | Eigenvalue Ai | Eigenvalue Az | Type of stability 
2 1/8000/4400;/- 3.49 42)]0 4542/Unstable 
2 2/8000/4600]- 3.59 68]0 43 6 8/Unstable 
2 3/8000/4800}]- 3.698 71)]0 418 7/Unstable 
2 4/8000;5000;/- 3.8000)]0 400 0/Unstable 
2 5}/8000/5200/-3.9007/)/0 3 80 7/Unstable 
2 6/8000}5400}-4.0007)]0 3 60 7;/Unstable 
2 7/8000;}5600}]- 4.100310 3 40 3}/Unstable 
2 8/8000}/5800}]- 4.199 3)]0 3 19 3/Unstabtle 
2 9/8000/6000}- 4.297 8]0 29 7 8}/Unstabtle 
3 0;/8000/6200]- 4.395 9)]0 275 9}/Unstable 
Table 8 Evaluating the impact of *2 between 6400 and 8200 on the type of stability 
Example | x 1 |X 2 | Eigenvalue Ai | Eigenvalue Az | Type of stability 
3 1/8000;6400|]-4.4935)/0 253 5/Unstable 
3 2/8000/6600/-4.5907]0 23 07;/Unstabtle 
3 3/8000/6800/-4.6875]0 207 5;/Unstabtle 
3 41/8000|}7000|-4.7839]/0 183 9/Unstable 
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3 5;/;8000]}7200);-4.8800]0 1600/Unstabltle 
3 6/8000)7400)|)-4.975710 13 5 7/Unstabtle 
3 7/8000)/7600)-5.0711 140 111%1/Unstabdtle 
3 8}/8000)/7800]-5.16611]0 008 6;Unstable 
3 9/8000)/8000j]-5.2608)1]0 060 8;Unstable 
4 0;/8000)8200)-5.3553 10 03 53;Unstabtle 


What contribution has the present analysis made over the earlier formulation? While the original formulation simply described 
the deterministic interaction between two competing species, the in-depth study of the impact of the co-existence steady state 
solution on the type of stability was not considered. In order to extend this idea, we have utilized a sound mathematical reasoning 
to investigate the effect of the co-existence steady state solution on the type of stability. On the basis of this present analysis, we 
can clearly mention that the present numerical analysis is a cutting edge contribution over the previous mathematically track-table 
and ecological track-table system of competing species. 


4. CONCLUSION AND FURTHER RESEARCH 


The key achievement of this important study depends on its significance contribution: we have utilized the method of numerical 
simulation to clearly study the effect of varying the co-existence of the steady state solution on the type of stability. In particular 
when the *1 population size is varied for a fixed *2 population size, we have found several significant dominant unstable steady 
state solutions that consequently change to a few instances of stable co-existence steady state solutions. This aspect of bifurcation 
analysis has the potential to influence ecosystem functioning, planning and sustainable development. The full details of other 
parameter variations in the context of pure competition and other types of interactions such as mutualism, commensalism and 
predation will be the subjects of our further research. 
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